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Abstract: We present a systematic study of generalized transverse- momentum depen- 
dent parton distributions (GTMDs). By taking specific limits or projections, these GTMDs 
yield various transverse-momentum dependent and generalized parton distributions, thus 
providing a unified framework to simultaneously model different observables. We present 
such simultaneous modeling by considering a light-cone wave function overlap representa- 
tion of the GTMDs. We construct the different quark-quark correlation functions from the 
3-quark Fock components within both the light-front constituent quark model as well as 
within the chiral quark-soliton model. We provide a comparison with available data and 
make predictions for different observables. 
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1 Introduction 

The investigation how the composite structure of a hadron, consisting of near massless con- 
stituents, results from the underlying quark-gluon dynamics, is a challenging problem, as 
it is of non-perturbative nature, and displays many facets. What are the longitudinal mo- 
mentum distributions of partons in a fast moving unpolarized or polarized hadron? What 
amount of transverse momentum do these partons carry, and how large is the resulting 
amount of orbital angular momentum? What is the spatial distribution of quarks inside a 
hadron as seen by a vector probe (coupling to the charge of the system) , by an axial vector 
probe (coupling to the axial charge), or even seen by a more complicated probe? 

The Generalized Parton Correlation Functions (GPCFs) provide a unified framework 
to address and quantify such questions. The GPCFs parametrize the fully unintegrated 
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off-diagonal quark-quark correlator, depending on the full 4-momentum k of the quark and 
on the 4-momentum A which is transferred by the probe to the hadron; for a classification 
see refs. [1, 2]. They have a direct connection with the Wigner distributions of the parton- 
hadron system [3-5], which represent the quantum mechanical analogues of the classical 
phase-space distributions. 

When integrating the GPCFs over the light-cone energy component of the quark mo- 
mentum one arrives at generalized transverse-momentum dependent parton distributions 
(GTMDs) which contain the most general one-body information of partons, corresponding 
to the full one-quark density matrix in momentum space. The GTMDs reduce to different 
parton distributions and form factors as is shown in figure 1. The different arrows in this 
figure represent particular projections in the hadron and quark momentum space, and give 
the links between the matrix elements of different reduced density matrices. 

Such matrix elements can in turn be parametrized in terms of generalized parton distri- 
butions (GPDs), transverse-momentum dependent parton distributions (TMDs) and gen- 
eralized form factors (FFs) . These are the quantities which enter the description of various 
exclusive (GPDs), semi- inclusive (TMDs), and inclusive (PDFs) deep inelastic scattering 
processes, or parameterize elastic scattering processes (FFs). At leading twist, there are 
sixteen complex GTMDs, which are defined in terms of the independent polarization states 
of quarks and hadron. In the forward limit A = 0, they reduce to eight TMDs which de- 
pend on the longitudinal momentum fraction x and transverse momentum k± of quarks, 
and therefore give access to the three-dimensional picture of the hadrons in momentum 
space. On the other hand, the integration over k± of the GTMDs leads to eight GPDs 
which are probability amplitudes related to the off-diagonal matrix elements of the parton 
density matrix in the longitudinal momentum space. After a Fourier transform of A_l to 
the impact-parameter space, they also provide a three-dimensional picture of the hadron 
in a mixed momentum-coordinate space [6-8]. The common limit of TMDs and GPDs is 
given by the standard parton distribution functions (PDFs) , related to the diagonal matrix 
elements of the longitudinal-momentum density matrix for different polarization states of 
quarks and hadron. The integration over x leads to a bilocal operator restricted to the 
plane transverse to the light-cone direction and brings to the lower plane of the box in Fig 1. 
The off-forward matrix elements of this operator can be parametrized in terms of so-called 
transverse-momentum dependent form factors (TMFFs). Starting from the TMFFs, we 
can follow the same path as in the case of the GTMDs, and at each vertex of the basis of 
the box of figure 1 we find the restricted version of the operator defining the distributions 
in the upper plane. Therefore, integrating out the dependence on the quark transverse 
momentum, we encounter matrix elements parametrized in terms of form factors (FFs), 
while the forward limit of TMFFs leads to transverse-momentum dependent spin densities 
(TMSD). Both FFs and TMSDs have the charges as common limit. 

Although a variety of models has been employed to explore separately the different 
observables related to GTMDs, a unifying formalism for modeling the GTMDs is still 
missing. In order to achieve that, we will exploit the language of light-cone wave functions 
(LCWFs), providing a representation of nucleon GTMDs which can be easily adopted in 
many model calculations. In order to simplify the derivation, we will focus on the three- 
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Figure 1. Representation of the projections of the GTMDs into parton distributions and form factors. 
The arrows correspond to different reductions in the hadron and quark momentum space: the solid (red) 
arrows give the forward Hmit in the hadron momentum, the dotted (black) arrows correspond to integrating 
over the quark transverse-momentum and the dashed (blue) arrows project out the longitudinal momentum 
of quarks. The different objects resulting from these links are explained in the text. 

quark (3Q) contribution to nucleon GTMDs, postponing to future works the inclusion of 
higher-Fock space components. In this way, we can express the GTMDs in a compact 
formula as overlap of LGWFs describing the quark content of the nucleon in the most 
general momentum and polarization states. Then, using the projections illustrated in 
figure 1, we can discuss the complementary aspects encoded in the different distributions 
and form factors. 

The plan of the paper is as follows. In section 2, we discuss the formal derivation of 
the LCWF overlap representation of the quark contribution to GTMDs, specializing the 
results to two light-cone quark models, namely the chiral quark-soliton model (xQSM) and 
the light-cone constituent quark model (LCCQM). In section 3, we focus the discussion on 
the TMDs, GPDs, PDFs, FFs and charges. In particular, we derive the general formulas 
obtained from the projections of GTMDs, and then we discuss and compare the predictions 
from both the xQSM and the LCCQM. In the last section, we draw our conclusions. 
Technical details and explanations about the derivation of the formulas are collected in 
three appendices. 

2 Formalism 

2.1 Parton Correlation Functions 

The maximum amount of information on the quark distributions inside the nucleon is 
contained in the fully- unintegrated quark-quark correlator W for a spin- 1/2 hadron [2-5], 
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defined as 

w[,{{P,k,A,N;rj) = \ j ^ e^^^ (p', A'|V^(-f )r W V'd A). (2.1) 

This correlator is a function of the initial and final hadron light-cone helicities A and A', 
the average hadron and quark four-momenta P = {p' + p)/2 and /c, respectively, and the 
four-momentum transfer to the hadron IS. = p' — p. In this paper, we choose to work in 
the symmetric light-cone frame, see figure 2. The corresponding kinematics is given in 
appendix A. The superscript F stands for any element of the basis {1, 75, 7^, 7'^75, io"'"^75} 




P-A/2 P + A/2 



Figure 2. Kinematics for the fully- unintegrated quark-quark correlator in a symmetric frame. 

in Dirac space. A Wilson line W = W(— |, ||n) ensures the color gauge invariance of the 
correlator, connecting the points — | and | via the intermediary points — | -|- 00 • n and 
I -|- 00 • n by straight lines. This induces a dependence of the Wilson line on the light-cone 
direction n. Since any rescaled four-vector an with some positive parameter a could be 
used to specify the Wilson line, the correlator actually only depends on the four-vector 

where M is the hadron mass. The parameter = sign(n°) gives the sign of the zeroth 
component of n, i.e. indicates whether the Wilson line is future-pointing [r] = +1) or 
past-pointing (r/ = —1). 

The quark-quark correlators defining TMDs, GPDs, PDFs, FFs and charges corre- 
spond to specific limits or projections of eq. (2.1). These correlators have in common the 
fact that the quark fields are taken at the same light-cone time = 0. Let us then focus 
our attention on the A;~-integrated version of eq. (2.1) 

wfiiP, X, k^,A, N;r])= f dk^ wf^iP, k, A, N; r?) 

(2.3) 



II 



(y, AX-|)F WV^d )|p, A)|^^^„, 



where we used for a generic four- vector a'^ = [a'^ ,a~ ,a±\ the light-cone components = 
(a" zb a'^)/\/2 and the transverse components a± = {a^,a?), and where x = k~^/P~^ and 
k_i_ are the average fraction of longitudinal momentum and average transverse momentum 
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of the quark. A complete parametrization of this object in terms of GTMDs has been 
achieved in ref. [2]. GTMDs can be considered as the mother distributions of GPDs and 
TMDs. Even though in the present paper we restrict our discussions to TMDs, GPDs, 
PDFs, FFs and charges, the formahsm can readily be applied to the case of GTMDs and 
will be the subject of an upcoming paper. 

2.2 Overlap Representation 

Following the lines of refs. [9, 10], we obtain in the light-cone gauge = an overlap rep- 
resentation for the correlator (2.3) at the twist-2 level. Moreover, this being a preliminary 
study, we restrict ourselves to the 3Q Fock sector^ (and therefore to the region < x < 1 
with ^ = — A+/2P+). Many theoretical approaches like e.g. [11] suggest that higher Fock 
components yield at best a 20% correction to quark observables with up and down flavors. 
We then write the correlator (2.3) as the following overlap 

,[r] r . . 1 



W^,^^{P,x,k±,A,N;r,) 



E /[dx]3[d2A;j3^(^)VV(OV'A/3(r-)M[W, (2.4) 



where the integration measures are defined as 

3 

[dx]3 



[d'A;^]3 




i=l 



(2.5) 



The indices (3 and f3' refer collectively to initial and final quark light-cone helicities Aj 
and A^, respectively. The tensor M^^^^'^ then represents the transition from the initial 
configuration to the final configuration of quark light-cone helicities and depends naturally 
on the Dirac operator F. The 3Q LCWF ipAjsii") depends on r which refers collectively to the 
momentum coordinates of the quarks in the hadron frame ki = {yi,Ki±), see appendix A. 
The function 5{k) selects the active quark average momentum 



5{k) = ^Q{x) 6{x - Xi) 5^2) (A?_j 



ki I 



i=l 



Since the labeling of quarks is arbitrary, we choose to label the active quark with i = 1 
and the spectator quarks with i = 2, 3. We can then write 



5{k) = 3 e(x) 5ix - xi) 5(2) (ki_ - k 



(2.6) 



and 



with mFI^'^ = uip',X')Tn{p,X) (2.7) 

2P+yr^ 

where u{p, A) is the free light-cone Dirac spinor. The light-cone helicity of a spectator 
quark is always conserved. 

^We consider here only transitions which are diagonal in both flavor and color spaces. Flavor and color 
indices are then omitted for clarity. 
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2.3 Helicity and Four- Component Bases 

There are only four twist-two Dirac structures rtwist-2 = {7^1 ■^f"'^'''75i ^o'^'''75) 7'''75}- They 
correspond to the four kmds of transition the hght-cone hehcity of the active quark can 
undergo, see e.g. [13-15] 

^[7+]A'A ^ ^A'A^ Mi'-'^^^']^'^ = {ajf\ mI^+^sI^'^ = (^73)^'^ (2.8) 

with (Tj the three Pauli matrices. For further convenience, we associate a four-component 
vector^ to every quantity with superscript F 

air] ^^u^ {a'',a\a^a^) = (a[^+U'^"'^^^Ut^'^'''^^U'^^^^l) • (2.9) 

With this notation, the correspondence (2.8) takes the simple form 

^uX'X ^ (^i')A'A^ (2.10) 

where a'^ = (l,o^). One can think of a'^^'^ = (a'^)^'^ as the matrix of a mere change of 
basis, the one being labeled by the couple AA' and the other by 

a"^ = Y^a''^'>^axx'. (2.11) 

A' A 

In the literature one often represents correlators in the helicity basis, i.e. in terms of 
helicity amplitudes 

VFa'A',AA= ^W^A'A^.AA', (2.12) 
where = gupcr'' with o"'^ = (1, —B). The symbols and satisfy the relations 

i a'^^'^a.^r' = S^:6^, ^Tr [a>^a,] = ^ J] a^"' V.aa' = (2.13) 

A' A 

We find however more convenient to work in the four-component basis. We then introduce 
tensor correlators 

Wf"" = ^Tr [af^W] = ^Yl ^''^'^'W^A'A- (2-14) 

A'A 

Helicity amplitudes and tensor correlators are related as follows 

W^A'A',AA = ^VF^"fTM'A^.AA', W^''' = ^ J] a''^^'^"^'^VFA'A',AA. (2.15) 

A'AA'A 

2.4 3Q LCWF from Constituent Quark Models 

So far, the exact 3Q LCWF ip\j3{r) cannot be derived directly from the QCD Lagrangian. 
Nevertheless, we can try to reproduce the gross features at low scales using constituent 
quark models. We focus here on two phenomenologically successful models which also have 
the advantage of incorporating consistently relativistic effects: the light-cone constituent 
quark model (LCCQM) [13-15] and the chiral quark-soliton model (xQSM) [16-24]. In the 

^Note this is not a Lorentz four-vector but Einstein's summation convention still applies. 
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Table 1. Expression of the 3Q LCWF in the LCCQM and the xQSM. In the LCCQM, m is the constituent 
quark mass, Mo is the free invariant mass of the 3Q state, and u} is the free quark energy. In the xQSM, 
= |'5±|, Mn is the soliton mass, and Ei^v is the energy of the one-quark discrete level. The functions 
f//(y, K±) and n±) are shown in appendix B. 



Model 






K± 




LCCQM 


i{j{r) 


m + yMo 




yMo - u 


xQSM 




f//iy,i^±) 




vMn - -Kiev 



LCCQM, one describes the baryon system in terms of the overlap of the baryon state with 
a state made of three free on-shell valence quarks. The 3Q state is however not on-shell 
M ^ TWq = X^jf^i) where uji is the energy of free quark i and M is the physical mass of 
the bound state. Since the exact baryon state is unknown, one approximates the overlap 
by a simple analytic function and fits the free parameters in order to reproduce at best 
some experimental observables, like e.g. the anomalous magnetic moment and the axial 
charge. In the xQSM quarks are not free but bound by a relativistic chiral mean field 
(semi-classical approximation). This chiral mean field creates a discrete level in the one- 
quark spectrum and distorts at the same time the Dirac sea. It has been shown that the 
distortion can be represented by additional quark-antiquark pairs in the baryon [20] . Even 
though the xQSM naturally incorporates higher Fock states, we restrict the present study 
to the 3Q sector. The inclusion of higher Fock states is postponed to a future work. 

Despite the apparent differences between the LCCQM and the xQSM, it turns out that 
the corresponding LCWFs are very similar in structure. In both models, the spin-fiavor 
part of the wave function is separated from the momentum part. Moreover, canonical spin 
and light-cone helicity are simply connected by an SU{2) rotation. We can take advantage 
of this similarity in structure and develop a formalism in terms of the generic 3Q LCWF 
(remember that (3 = {Aj}) 

3 

^|;J,p{r) = Af^{r)Y,^'^r''''ll^x,'^^^h), (2.16) 

o"i i=l 

where M is a (real) normalization factor, ^(r) is a symmetric momentum wave function, 
(j)CTio-20-3 -g ^j^^ SU{6) spin-flavor wave function, and D{k) is an SU{2) matrix relating 
light-cone helicity A to canonical spin a 




Kn,L = K^±iKy. (2.17) 

In order to specify the forthcoming expressions to one of the models, one has to perform 
the substitutions given in table 1. In the LCCQM, the quarks are free and the matrix D{k) 
then simply corresponds to the Melosh rotation R^f [25] 

D\a[k) := DJ)^ {Rcf{k)) = , = . (2.18) 

' {m + yMo)'^ + 
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In the xQSM, the quarks are bound and the matrix D(k) is related to the one-quark 
discrete-level wave function F\^{k) created by the mean field [21] 



\K\D^,{k):=F^,(k). (2.19) 

The relation between light-cone helicity and canonical spin then involves the dynamics 
which is encoded in the functions f//{y,n±) and f±{y,K±). Note that in both models the 
rotation is about the axis^ k± x e^. For more details on the 3Q LCWF in the LCCQM 
and the xQSM, see appendix B. 

Connecting light-cone helicity to canonical spin (and consequently light-cone wave 
functions to instant-form wave functions) is usually an extremely difficult task since it 
involves boosts which contain the interaction. Only in simplified pictures is this connection 
tractable. In the two models we consider here, quarks are free or interact with a relativistic 
mean field, i.e. they are quasi- independent. The main effect of boosts is creating an angle 
between instant-form and light-cone polarizations. This angle vanishes when the particle 
has no transverse momentum. An equivalent point of view is to say that a quark state 
with definite light-cone helicity corresponds to a linear combination of quark states with 
both canonical spin '[ and J,, responsible for the non-diagonal elements in D(k). 

2.5 3Q Proton Amplitude 

Inserting the LCWF given by eq. (2.16) in the overlap representation of the correlator 
tensor W^^'^ (2.4), we obtain 

where A^'^{r',r) stands for 

^^-(r', r)=A Of ih ■h) + B [/^ {h ■ O^f + [h ■ OiY] . (2.21) 

The coefficients A and B are flavor factors. We used = O!^^ with the matrix O^*^ given 
by 

/ K' ■ K i (^K' X k\ i( K' X k\ -i (k' X k\ \ 



\K'\\K\ 



iiK' X K] K' -K- 2K'^K^ -K^y - K'yK^ K'^K^ + K'^K^ 



[k' X -K'yK^ - K'^Ky K' -K- 2K'yKy K'yK, + K'.Ky 

y (^K' X -K',K^ - K'^K^ -KKy - K'yK^ -K' ■ K + 2K'^K^j 

(2.22) 

Details on the derivation can be found in appendix C. 

Let us interpret this master formula. The correlator tensor W^^ has two indices ^ 
and V which refer to the transition undergone by the hadron and active quark light-cone 



This is not surprising by noticing that the generators of transverse boosts on the light cone are given 
by — + J1 X , where B± and Jl are the ordinary boost and rotation generators [26]. 
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helicities, respectively, eq. (2.20) expresses the tensor correlator in terms of the 3Q overlap 
/[da^]3 [d^/c_L]3 of initial ^(r) and final ^*(r') symmetric momentum wave functions 

with the tensor A^^{r' ^r^ for a fixed mean momentum of the active quark A(A;). The 
tensor A^"{r' corresponds to the overlap of the three initial quarks with the three final 
ones. In the S'?7(6) spin-flavor wave function cj)^i'^2 0-3 quark canonical spins is 

always aligned with the hadron helicity A. Hence A^^ is the sum of two contributions: the 
hadron helicity can be aligned with the canonical spin of either the active quark i = 1 or 
one of the spectator quarks i = 2, 3. Since the spectator quarks are equivalent, they enter 
in a symmetric way in eq. (2.21). The coefficients A and B then give the weight of each 
contribution and depend on the flavor of the active quark and the nature of the spin-1/2 
hadron. In a proton, we have 

Al = ^, B; = 1, ^^ = -1, S^ = 2. (2.23) 

Finally, the matrix O^'^ in eq. (2.22) describes the overlap of an initial quark with a final 
quark. Rows and columns correspond to the type of transition undergone by the quark 
canonical spin and light-cone helicity, respectively. 

For example, consider the vector operator F = 7+. eq. (2.8) tells us that this operator 
is not sensitive to the active quark light-cone helicity. To see what happens in terms of 
active quark canonical spin, we have to look at the first column of eq. (2.22). It appears that 
the vector operator is generally sensitive to the canonical spin {O^^ 7^ for j = 1,2,3). 
Note that in absence of momentum transfer K' = K, the vector operator becomes also 
insensitive to the active quark canonical spin. This can be easily understood as follows. 
The orientation of canonical polarization relative to light-cone polarization depends on the 
momentum of the particle. In presence of momentum transfer, the initial and final SU (2) 
rotations D{k) are different and the vector operator effectively sees the difference between 
canonical spin t and |. 



3 Results and Discussion 

We now apply and specialize the general formalism described in the previous sections to 
the extraction of TMDs, GPDs, PDFs, FFs and charges. The normalization constant M 
for each model has been fixed so as to satisfy the valence-quark sum rules. Namely, there 
are (in an effective way) two up quarks and one down quark in a proton. In the xQSM, 
the one-quark discrete level wave function is given by the sum of a bare discrete level 
contribution F^^^ and a relativistic contribution F^^^ due to the distortion of the Dirac sea. 
The latter one is usually not discussed because it is expected to give only a small correction. 
Here we calculate this contribution for the first time, as explained in appendix B, and we 
explcitly check that it is in general not relevant. Therefore, except for the charges, all the 
results in the xQSM will refer only to the F^^^ contribution. 
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3.1 Transverse Momentum-Dependent Distributions 



The forward limit A = of the correlator W is given by the quark-quark correlator, denoted 
as ^> 



cI>P^(P, X, k^, N; r?) = W^^iiP, x, , 0, N; r?) 

1 fdz-d'^z± 



(27r)3 



ixP^z —ik±-z± 



(p,A'|V'(-i)rwv(f)|P,A) 



^=0 

(3.1) 



It is parametrized by TMDs at leading twist in the following way 



M "-1 



ky ^_|_ 

M JIT 
_ kx f± 
M JIT 





^1 + ^aF 



IT 



TV?2- '"IT 



kx l,± 

M 

-mt n 

~2M^ 



IT 



IT 



«■!/ /,-L 



\ 

TtaiT 

9lL / 



(3.2) 



TMDs are functions of x and k'j_ only. The multipole pattern in k± is clearly visible in 
eq. (3.2). The TMDs /i, giL and /ii give the strength of monopole contributions and 
correspond to matrix elements without a net change of helicity between the initial and 
final states. The TMDs fy^p, qit, hj^, and hj^^ give the strength of dipole contributions 
and correspond to matrix elements involving one unit of helicity flip, either on the nucleon 
side {fyj^ and gix) or on the quark side {hj^ and hj^^). Finally, the TMD hj^p gives the 
strength of the quadrupole contribution and corresponds to matrix elements where both 
the nucleon and quark helicities flip, but in opposite directions. Conservation of total 
angular momentum tells us that helicity flip is compensated by a change of orbital angular 
momentum [27] which manifests itself by powers of k±/M with M the mass of the nucleon. 
In the spirit of [28], we define the transverse (0)-, (1/2)- and (l)-moments of a generic 
TMD j{x,kj_) as 

j(0)(x)^ I d^k^j{x,kl), 
3^"'\^)^ ld'k^^^j{x,kl), (3.3) 
J«(x) = ld\^^^j{x,kl). 



The transverse (1/2)- and (l)-moments are chosen such that they directly represent the 
strength of the dipole and quadrupole distributions as function of x, respectively. Note 
that our definition of the (l/2)-moment is twice larger than in [28]. 
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In absence of momentum transfer, i.e. A = 0, the matrix O'^'^ in (2.22) reduces to^ 



Qf,^ A=o J_ 
i?2 








_ -2K^Ky 



2K^K, 



2KyK, K^-2Kl 2KyK, 
\ -2K,K, -2K,Ky -i?2 + 2KlJ 



(3.4) 



and the tensor A^'^ of eq. (2.21) becomes 

^00 = + 2B) , A^^ = = 0, 



(3.5) 



Note that A^^ = = because we are not considering gluon degrees of freedom, i.e. the 
generic quark wave function (2.16) leads to vanishing Sivers and Boer-Mulders functions. 
Moreover, the distributions for different flavors are just proportional, which is a consequence 
of the underlying SU{&) spin-flavor symmetry. Let us then introduce the spin-flavor factors 
A^'' and P"^ for a flavor q as follows 



m = {A'' + 2B'') /3 and = Ay3. (3.6) 

With this definition N'^ and P"^ can be identified, respectively, with the number of quarks 
of flavor q in the baryon and the non-relativistic contribution to the total spin of the baryon 
coming from quarks of flavor q. For a proton, we have 

N^ = 2, A^p=l, Pp" = 4/3, Pp^ = -l/3. (3.7) 

Our model results for some transverse moments of TMDs are shown in Fig 3. The 
shape of the curves within the xQSM and the LCCQM are very similar. The size of 
the longitudinal and transversity distributions is somewhat smaller for the LCCQM. On 
the other hand, the peak of the distributions for the transverse moments of the polarized 
TMDs is larger in the LCCQM than in the xQSM, especially for hf^\ This pattern for 
the magnitude of TMDs in the two models indicates that there is more orbital angular 
momentum in the LCWF of the LCCQM than in the xQSM. This can be understood as 
follows. In these two models, the instant-form wave functions do not contain any orbital 
angular momentum. The proton spin originates solely from the quark spins. However, the 
corresponding light-cone wave functions do involve orbital angular momentum since quark 
canonical spin and light-cone helicity do not coincide in general. The proton spin originates 
from both the quark light-cone helicities and orbital angular momentum. In other words, 
all the orbital angular momentum that appears in this paper comes from the generalized 
Melosh rotation of eq. (2.17), or more precisely from its non-diagonal elements K±. As 
one can see from eqs. (3.2) and (3.4), all the functions {giL — hi), qit, hj^^ and hj?p vanish 
identically in the limit K± — )• 0, i.e. in absence of orbital angular momentum. It is then 

^This matrix is orthogonal and composed of two blocks O = ( J jj ) where R is an 50(3) matrix. This is 
hardly surprising since in this case the transformation O — {k)M D{k) is just the well-known homomor- 
phism between SU{2) and SO{3). 
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Figure 3. Some examples of transverse moments of TMDs as function of x (see text for their definition). 
In all panels the results are for the "flavorless" TMDs. TMDs of definite flavor follow from multiplication 
by the spin-flavor factor A'^'' in the unpolarized case and in the polarized case (see eq. 3.7). Solid curves: 
results in the ;;j;QSM. Dashed curves: results from the LCCQM of ref. [27]. All the results are at the 
hadronic scale of the models. 



hardly surprising that the TMDs we obtained are not ah independent. There exist three 
relations among polarized TMDs, which are flavor independent and read 



giT + Kl = 0, (3. 



9lL 



0, (3.9) 



giT + 2hihtT = 0. (3.10) 



One further flavor-dependent relation involves both polarized and unpolarized TMDs, and 
is given by 

PVi' + 5?i = 2/i?, (3.11) 

where = P'^/N'^. For example, from the relation (3.9) we see that the difference between 
helicity and transversity distributions qil — hi is directly connected to the pretzelosity 
distribution /i^. The difference being smaller in the xQSM, it follows immediately that 
the pretzelosity distribution is also smaller compared to the LCCQM. 

The relations in eqs. (3.8)-(3.11) actually hold in a large class of quark models (see 
refs. [28, 29] and references therein). Their physical origin has been briefly discussed in 
[30] and will be explained in more details in a forthcoming publication. Although such 
quark model relations are appealing, one should keep in mind that they break down in 
models with gauge-fleld degrees of freedom and are not preserved under QCD evolution. 
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Despite these limitations, they can provide useful guidelines for TMD parametrizations 
to be further tested in experiments. An interesting result comes from recent lattice cal- 
culations [31, 32], which give hf^^ ~ ~9rT' favor of the relation (3.8). In particular, 
they obtained for the average dipole deformation (along the direction of the nucleon po- 
larization) of the k± density for longitudinally polarized quarks in a transversely polarized 
nucleon {k^igix)) = M^i^^"//^"^" = 67(5) MeV, for up quarks, and (k^igir)) = -30(5) 
MeV, for down quarks. The corresponding values for transversely polarized quarks in a 
longitudinally polarized nucleon are {k^ihf-i)) = M hfj^^^ ^ / f[^^ ^ = —60(5) MeV, for up 
quarks, and {k^{hj^j^)) = 16(5) MeV, for down quarks. These results are remarkably sim- 
ilar to our quark model calculations: {k^{giT)) = — (^"(^li)) ~ 54.10 and 55.8 MeV, 
and {kiigir)) = -{kHh^L)) = -27.05 and -27.9 MeV in the xQSM and the LCCQM, 
respectively. 

Note also that, in the xQSM, the structure functions do not vanish at x = while they 
do vanish in the LCCQM. In the LCCQM, the functional form of 3Q LCWF is assumed on 
the basis of phenomenological arguments. Using the simple power-law Ansatz, the LCWF 
itself vanishes when any yj — )■ (see appendix B). This is at variance with the wave 
function of the xQSM which comes from the solution of the Dirac equation describing the 
motion of the quarks in the solitonic pion field. We should also mention that the xQSM 
was recently applied in ref. [33] to calculate the unpolarized TMDs, taking into account 
the whole contribution from quarks and antiquarks without expansion in the different Fock 
components, and obtaining results for the isoscalar combination u + d which corresponds 
to the leading contribution in the l/Nc expansion. 



3.2 Generalized Parton Distributions 

When integrating the correlator W over fcj^, one obtains the quark-quark correlator, de- 
noted as F 



f[,{{P, X, A,N)= I d^k^ W[,{{P, X, k^,A, N; ??) 



1 f dz 
2 



AxP+z 1^1 \l 



2-K 



(p',AX-|)rwV'(i)|p,A) 



(3.12) 



Z+=2j^=0 



Note that the integration over kx_ removes the dependence on r/, and we are left with 
a Wilson line connecting directly the points — | and | by a straight line. Furthermore, 
working in the light-cone gauge = 0, the gauge link can be ignored. The correlator F 
is parametrized by GPDs at leading twist in the following way 
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(3.13) 
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where we used the "natural" combinations of standard GPDs^ 



e 



Et + 



Ht + C Et 



1-e 



£ 
£ 

£t 
£t 



E 



iE 



IHt + Et — ^ Et 
Et — ^ Et 



(3.14) 

Comparing eq. (3.13) with eq. (3.2), one notices a strong analogy. Essentially, the same 
multipole pattern appears, where the role of k±_ in eq. (3.2) is played by A_l in eq. (3.13). 
Going to impact-parameter space representation does not change the multipole structure 
[34]. This suggests that there might be relations or dynamical connections between GPDs 
and TMDs, e.g. between E and or between 2Ht + Et and hj^ [35, 36]. There is however 
no direct link since GPDs and TMDs often originate from different mother distributions 
[2, 37], but one might still expect some correlation between signs or similar orders of 
magnitude from dynamical origin (there are after all deep connections with quark orbital 
angular momentum). 

Results for the GPDs as function of x and different values of and t are shown in 
figures 4-7. Since we are considering only quark degrees of freedom, the calculations are 
restricted to the DGLAP region for x > ^. The xQSM was already applied to obtain 
predictions for the GPDs within a different framework [38-44], i.e. using the instant-form 
quantization and incorporating the full contribution from the discrete level and Dirac sea, 
without expansion in the different Fock-space components. In this way the full range 
of X was explored, but only for the flavor combinations which correspond to the leading 
contribution in the 1/Nc expansion. 

For all the GPDs we show the separate up and down-quark contributions. Note that 
in contrast to TMDs, GPDs for different quark flavors are not simply proportional. As 
long as the momentum transfer t = is not vanishing, the two terms in eq. (2.21) are 
usually different and non-zero. 

The point x = ^ corresponds to vanishing longitudinal momentum for the active quark 
in the final state (see eq. (A. 2) in appendix A). Therefore, similarly to the case of the TMDs, 
it is a node for the wave function in the LCCQM, but not for the xQSM. Accordingly, all 
the GPDs from the LCCQM are vanishing in this point. 

The distributions of unpolarized quarks are described by the GPDs H and E, the last 
one being non-vanishing only in the presence of orbital angular momentum. Comparing 
the results from the LCCQM and the xQSM in figure 4, we see that they are very similar 



^The fundamental physical object is the matrix element. The GPDs are defined according to a specific, 
but not unique, parametrization of matrix elements. Therefore they have no a priori simple interpretation. 
Combinations at ^ = appeared already in [34] . 
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Figure 4. Results for the spin averaged {H'', three upper panels) and the helicity flip three lower 
panels) generalized parton distributions for the up (left panels) and down (right panels) flavors, at fixed 
values of ^ and t as indicated. Solid curves: results in the xQSM. Dashed curves: results from the LCCQM 
of ref. [13]. 

in size and in the behavior at large x. The faster fah off of with respect to H'^ for x — t- 1 
is due to the decreasing role of the orbital angular momentum for increasing longitudinal 
momentum of the active quark. This feature is common also to the other GPDs which 
arise from a transfer of orbital angular momentum between the initial and final state, like 
the chiral-odd GPDs Et and Et shown in figures 6 and 7, respectively. For longitudinally 
polarized quarks, we only show the results for the GPD in figure 5. We refrain from 
presenting the quark contribution to the E"^ without discussing also the pion-pole term. 
As a matter of fact, the pion-pole is by far the largest contribution to this GPD and the 
differences between the LCCQM and the xQSM for the quark contribution would not be 
relevant for the total results. For the chiral-odd GPDs in figures 6 and 7, at ^ = there 
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Figure 5. Results for the helicity-dependent generalized parton distribution H'^ with the same notation 
as in figure 4. Results for the LCCQM are from ref. [14]. 



is no Ej, because it vanishes identically being an odd function of ^ as consequence of 
time-reversal invariance. 

Another common property to the model results for all GPDs concerns the t-dependence. 
In general it is more pronounced in the low x region and it affects the position of the peak, 
especially in the LCCQM calculation where we observe a shift to larger x for increas- 
ing values of t. On the other hand, the distributions show a very weak i-dependence at 
large x < 1, and are all vanishing at the end point j; = 1 as expected from momentum 
conservation. 



3.3 Parton Distribution Functions 

Considering the forward limit A = of the correlator F or, equivalently, integrating the 
correlator $ over the quark transverse momentum k± yields the quark-quark correlator 



Tl^^{P,x,N) = Fj^,[{P,x,0,N) 

d^k^wl}jP,x,k^,0,N;7]) 



(3.15) 



2./ 27r 
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z+=zj_=0 
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Figure 6. Results for the ehiral-odd H!^ (three upper panels) and E!^ (three lower panels) gener- 
alized parton distributions with the same notation as in figure 4. Results for the LCCQM are from 
ref. [15]. 



It is parametrized by PDFs at leading twist in the following way 



//i 0\ 

/ii 

/ii 

\0 51/ 



(3.16) 
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Figure 7. Results for the chiral-odd H!^ (three upper panels) and E!^ (lower panels) generalized 
parton distributions with the same notation as in figure 4. Results for the LCCQM are from ref. [15]. 



Naturally, only the monopoles in k± and A_l survive and one obtains the well known 
relations between PDFs, GPDs and TMDs 

h{x) = H{x,0,0) = I d^k^fi{x,kl), 

gi{x) = Hix,0,0) = J d?k^giL{x,kl), (3.17) 
hi{x) = Ht{x,^,Q) = y hi{x,ki). 

In figures 8-10 we show the results for the unpolarized, helicity and transversity dis- 
tributions from the xQSM and the LCCQM after appropriate evolution from the hadronic 
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Figure 8. Results for the unpolarized PDF xfi for up (left panel) and down (right panel) quark. The 
solid (dashed) curves are the results from the xQSM (LCCQM) after NLO evolution from the model scale 
Qo = 0.259 GeV^ to = 5 GeV^ The crosses are the fit to the experimental data from the CTEQ 
analysis at NLO of ref. [45]. 




X X 



Figure 9. Results for the polarized PDF xgi for 
solid (dashed) curves are the results from the xQSM 
Ql = 0.259 GeV^ to = 5 GeV^ The crosses are 
Ref. [46]. 



up (left panel) and down (right panel) quark. The 
(LCCQM) after NLO evolution from the model scale 
the fit to the experimental data from the analysis of 



scale of the models to the relevant experimental scales. A key question emerging not 
only here but in any nonperturbative calculation concerns the scale at which the model 
results for the parton distributions hold. From the point of view of QCD where both 
quarks and gluon degrees of freedom contribute, the role of low-energy quark models is 
to provide initial conditions for the QCD evolution equations. Therefore, we assume the 
existence of a low scale Qq where glue and sea-quark contributions are suppressed, and 
the dynamics inside the nucleon is described in terms of three valence quarks confined by 
an effective long-range interaction. The actual value of Qq is fixed by evolving back the 
unpolarized data, until the valence distribution matches the condition that its first moment 

< X >y is equal to the momentum fraction carried by the valence quarks as computed in 
the model. Since in our models only valence quarks contribute, the matching condition is 

< x{Qq) >y= 1. Starting from the initial value < x{Q'^) >t,w 0.36 at = 10 GeV^ and 
fixing the values of Aqcd and heavy-quark masses as in Ref. [45], we find QqIlo = 0.172 
GeV^ and QqInlo = 0.259 GeV^ after LO and NLO backward evolution, respectively [47]. 
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Figure 10. Results for the transversity distribution xhi for up (left panel) and down (right panel) quark. 
The solid (dashed) curves are the results from the xQSM (LCCQM) after NLO evolution from the model 
scale Qo = 0.259 GeV^ to — 2.5 GeV'^. The shaded area corresponds the uncertainty band due to the 
statistical error of the parametrization of refs. [48, 49]. 



In the case of the unpolarized and polarized PDFs in figures 8 and 9, respectively, we show 
the results only for the non-singlet (valence) contribution after NLO evolution to = 5 
GeV^, since the models at the hadronic scale consider valence quarks, and the gluon and sea 
contribution at higher scales are generated only perturbatively. The results of the xQSM 
and the LCCQM are very similar, for both up and down quarks, and overall reproduce the 
behavior of the PDFs as obtained from phenomenological parametrizations for unpolar- 
ized [45] and polarized [46] PDFs fitted to experimental data. However, we notice a faster 
falloff of the tail of the up-quark distributions at larger x in our models with respect to the 
parametrizations. For the transversity distribution in figure 10, we show the model results 
after NLO evolution to = 2.5 GeV^, in comparison with the available parametrization 
from refs. [48, 49]. In this case, since there is no gluon counterpart and the sea can not be 
generated perturbatively, the model results after evolution correspond to the pure valence 
contribution. Compared with the phenomenological parametrization, our results are larger 
for both up and down quarks. However, we find that hi is smaller than the Soffer bound [50] 
calculated from the model predictions for fi and gi, i.e. < ^[fi{x) + gl{x)]. The 

results within the LCCQM for the transversity were also used in ref. [51] to predict the 
Collins asymmetry in semi-inclusive deep inelastic (SIDIS) scattering. By using the Collins 
fragmentation function H^- of ref. [54], the work of ref. [51] obtained a very good agreement 
with available HERMES [52] and COMPASS [53] data. However other extractions of 
in the literature [48, 49, 54, 55] are based on different assumptions. In particular, issues 
concerning evolution effects [56] in the extraction of the Collins function from data span- 
ning a large range of are not yet settled. The Collins function from [54] was extracted 
from a fit to SIDIS data [52, 53]. On the other hand, the authors of Refs. [48, 49] per- 
formed a simulatenous fit of SIDIS and e~^e~ annihilation data from BELLE [57], including 
approximately effects of the evolution in Q^. The difference in the results for from the 
two analyses explains why both the results for hi in figure 10, within the LCCQM and 
from the extraction [48, 49], are consistent with the data. 

We also mention that the calculation of the helicity and tansversity distribution in the 
xQSM, without expansion in the Fock space and using instant-forn quantization, can be 
found in references [58, 59]. 
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3.4 Form Factors 

When integrating the correlator F over x, one obtains the quark-quark correlator A 



1 



2P^ 



A'AV 

(p',A'|v^(o)rv(o)|p,A). 



(3.18) 



Since in this case we deal with a local operator, the Wilson line drops out together with 
the dependence on the light-cone direction n. This correlator A is parametrized by FFs in 
the following way 
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(3.19) 



FFs appear as the lowest x-moment of GPDs and are scale independent 
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(3.20) 



They are also independent of ^ as a consequence of Lorentz invariance (polynomiality of 
GPDs). There is no FF associated to Et because time-reversal invariance implies that Et 
is odd in ^, i.e. E^ix, —S,,t) = —ET{x,S,,t). By polynomiality, its lowest x-moment must 
then vanish^ [60]. For illustration, the calculated nucleon electromagnetic Sachs FFs 



GEiQ^) = Fi(Q') - tF2(Q2), Gm(Q') = Fi(Q') + F2iQ^) 



(3.21) 



with T = Q'^ /AM"^ are compared with existing experimental data in figure 11. In particular 
we plot the results for Gm/{^Gd) of the proton and neutron (upper panels), with the 
standard dipole form factor Gr, = 1/{Q'^ + A|,)^ and A^, = 0.71 GeV^, and the results for 
pPG^^/G^j and G^ (lower panels), for the proton and neutron, respectively. The values for 
the magnetic moments are those of the models (see table 2 and subsequent discussion). 



®Note that the matrix element from which Gp can be extracted has an explicit ^ factor. By analogy, 
we might define formally a fourth tensor structure function as H4{t,^) — J^-^ dx ET{x,^,t)/£, which is not 
forced to vanish. It does however not correspond to any known Lorentz structure in the parametrization 
of the correlator A. In this case, the polynomiality argument does not apply a priori, and H4 might also 
depend on 5. 
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Figure 11. The four nucleon electromagnetic Sachs FFs compared to the world data: MIT (purple stars), 
MAMI (red circles), JLab (blue squares), as well as older results (black triangles and open symbols). The 
references to the data can be found in [61]. Solid curves: results in the xQSM. Dashed curves: results from 
the LCCQM of Ref. [62]. The values for the anomalous magnetic moments are given in table 2. 



The LCCQM reproduces rather well the trend of the proton data up to w 1 
GeV^. At higher values of the momentum transfer, the slope of is too steep and 
the curve deviates from the data. This effect is somehow compensated if we look at the 
ratio /U^G^/G^. Here the combined effect of a slightly overestimated G^j and a slightly 
underestimated G^ gives a result for the ratio which follows the fall-off of the data up to 
~ 2 GeV^. In the neutron case, the description of the form factors within the LCCQJVI 
is less satisfactory. In particular G^ is largely underestimated and we are not able to 
reproduce the slope at low Q^. The neutron radius originates as a partial cancellation 
between the Dirac radius and the contribution of the anomalous magnetic moment (the 
so-called Foldy term) which dominates. In the LCCQM, the Dirac radius is much larger 
as compared with the xQSM, leading to a larger cancellation in the neutron radius, which 
reaches only 2/3 of the value in the xQSM. As it was shown in Ref. [62], the description in 
the LCCQM can be improved by taking into account the contribution of the meson cloud 
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of the nucleoli and by relaxing the approximation of SU (6) symmetry in the model. The 
combined effects of a small percentage of mixed-symmetric terms in the nucleon LCWF 
and of the pion-cloud, which mainly acts at very small Q^, provides a much improved 
understanding of within that model. Those effects are much less significant on the 
other nucleon form factors. 

The results within the light-cone xQSM exhibit a peculiar behavior for the magnetic 
form factors. Both in the proton and neutron case the results rise faster than the dipole 
form at low Q^, while have a steeper fall off at higher Q^. The deviation is within 10%, but 
it goes in the opposite directions in the two ranges of < < 1 GeV^ and 1 < < 2 
GeV^, in such a way that the global curvature is quite different from the data. On the 
other hand, the results for are in much better agreement with the data. The slope at 
low is much steeper than in the LCCQM, going in the direction of the experimental 
data which support even a steeper rise. Finally, the results for ^^G^^/G^^ have a too fast 
falloff with Q^, and deviate from the experimental data by ~ 20% in the whole range of 

3.5 Charges 

Considering the forward limit A = of the correlator A or, equivalently, integrating the 
correlator IF over x yields the quark-quark correlator Q 

dxFl}^{P,x,N) 



2P 

which is parametrized by the charges in the following way 
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(3.23) 



The vector charge q = Ff{0) = J dx fKx) represents the (effective) number of quarks 
with flavor q. The normalization constant M of the models has been chosen such that 
u = 2 and d = 1 for the proton. The axial charge Aq = G^(0) = J dx gKx) represents the 
fraction of baryon helicity carried by the spin of quarks with flavor q. Finally, the tensor 
charge dq = Hf(0) = J dxh1{x) represents the fraction of baryon transversity due to the 
spin of quarks with flavor q. From eq. (3.4), one can see that the axial and tensor charges 
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Table 2. Results for the axial charge Aq, tensor charge 5q, anomalous magnetic moment n'', and 
tensor anomalous magnetic moment k'^ in the xQSM and the LCCQM, compared with experimental and 
phenomenological values. xQSM I refers to the calculation with F^^^ only, while the results labeled xQSM II 
and III include the effects of the Dirac sea on the discrete level, without and with Pauli-Villars regularization, 
respectively (see appendix B). The results within the LCCQM are from refs. [62-64]. The experimental 
data for the axial charges at = 5 GeV^ are from [65] and the anomalous magnetic moments are from [66]. 
The model results for the tensor charges are evolved at LO to — 0.8 GeV^ for comparison with the 
values from the phenomenological extraction of Ref. [49] . 
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-0.444 


*J.04„Q 22 


n 00+0.09 
"•^"-•-0.16 


1.673 


-2.033 







are usually different. Only in the non-relativistic limit, i.e. when there is no appreciable 
orbital angular momentum kj^ w and thus no distinction between canonical spin and 
light-cone helicity, they do coincide. This is the origin of the claim that the difference 
between axial and tensor charges is a measure of quark orbital angular momentum. One 
should however keep in mind that this is not completely true in general because of higher 
Fock components. The FFs F2{t), H2{t), H3{t) and Gp{t) being associated with dipole 
or quadrupole structures in A_|_ do not have corresponding charges. Nevertheless, one can 
still define their limit for vanishing momentum transfer t = (as long as it is finite). For 
example, the anomalous magnetic moments and anomalous tensor magnetic moments are 
defined by k'^ = -F|(0) and = H2{0), respectively. 

In table 2, we compare the calculated charges and anomalous moments for a proton 
with the corresponding experimental values and phenomenological extractions. For the 
xQSM we quote the results with the bare discrete-level contribution only (xQSM I) and 
with the relativistic contribution coming the distortion of the Dirac sea, without (xQSM 
II) and with Pauli-Villars regularization (xQSM III), see appendix B for more details. The 
corrections due to the effects of the sea are in general small, especially after regularization. 
Therefore, in the following discussion of the results we will refer to the values of xQSM I 
only. 

Since the isovector combination of the axial charge qa = Att — Ad is scale independent 
(as a consequence of current conservation) and since the isoscalar combination Au + Ad 
is weakly scale dependent, we did not evolve the model results for the axial charges. The 
last ones and the anomalous magnetic moments tend to overestimate the contribution from 
up quark and, vice-versa, to underestimate the down quark contribution. The deviations 
are such that they compensate in the results for the isovector axial charge g^- We find 
gA = 1-24 and gA = 1-44 for the LCCQM and the xQSM, respectively, to be compared with 
the experimental value gA = 1-269 obtained from the HERMES data [65], in agreement 
with the value obtained from /3-decay measurements [66]. On the other hand, for the 
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isoscalar contribution Au + Ad, the deviations go in the opposite direction, such that 
about 75% (85%) of the nucleon hehcity is carried by the quark hehcities in the LCCQM 
(xQSM). This has to be contrasted with the much smaller experimental result of about 
37% at a scale around 10 GeV^ [67]. 

For the anomalous magnetic moments, the deviation is bigger for down quarks than 
for up quarks. This mainly affects the neutron results, for which the discrepancy from the 
experimental data is bigger than in the proton case. In particular, we find /i^ = 2.78 and 
/x" = -1.69, within the LCCQM, and = 2.69 and /i" = -1.62, within the xQSM, to 
be compared with the experimental values /.t^ = 2.793 and /.t" = —1.913. A significant 
improvement is obtained for in the LCCQM by taking into account the meson-cloud 
contribution as in Ref. [62]. This is at the price of slightly worse results for gA_ and 
fiP. As a matter of fact, it is in general a challenging task in light-cone quark models to 
reproduce simultaneously all the three quantities, see e.g. refs. [68-71]. While it has already 
been shown in refs. [21-23] that higher Fock components within the xQSM contribute 
significantly to the axial charges, it remains to be checked explicitly that the same holds 
for the magnetic moments. 

For the tensor charges, we evolved the model results at LO to = 0.8 GeV^ for 
comparison with the corresponding results of the phenomenological extraction of Ref. [49]. 
In this case, both the LCCQM and the xQSM agree very well with the phenomenological 
value for the down quark, while the up-quark contribution is larger by ~ 40%. We note 
however that the values of the quark tensor charges are strongly scale dependent and may 
depend crucially on the choice of the initial scale of the models [42]. A safer quantity to 
compare with is the ratio 6d/5u which is scale independent. In our models, the assumption 
of SU{6) symmetry implies 5d/6u = —1/4, which is compatible, within error bars, with 
5d/5u = — 0.42;[Iq2o'^^ from the fit of [49]. We note that the deviation of the experimental 
data from the SU{6) limit of 1/4 for the ratio Ad/Au is much more significant. The 
experimental data at = 5 GeV^ give Ad/Au = —0.54, and this value is weakly scale 
dependent. Therefore we expect that the inclusion of SU{6) symmetry breaking terms can 
play an important role for the description of the longitudinal spin degrees of freedom, while 
it is less relevant for observables related to the transverse spin. Note that in the xQSM the 
SU{6) symmetry holds only at the 3Q level. Let us mention that the instant-form version 
of the xQSM of Ref. [72], which includes in principle all Fock components, gives the ratios 
5d/du = -0.3 and Ad/Au = -0.45 at the model scale Ql = 0.36 GeV^. 

Finally, the results for the tensor anomalous magnetic moments are very similar in 
the LCCQM and the xQSM. The values of k^/2M and k^/M are a measurement of the 
average distortion of the density describing the distribution in impact-parameter space 
of transversely polarized quarks in unpolarized nucleons. Therefore, although the model 
results give Kj. > Kj^, the distortion in the spin density is larger for down quarks than for 
up quarks by about 30%. The corresponding quantity in the transverse- momentum space 
is given by the average dipole distortion induced by the Boer-Mulders function. In this 
case, using the recent LCCQM calculation of Ref. [73], we find the same results for the 
relative size of the distortion in the down-quark distribution with respect to the up-quark 
distribution. The same arguments apply for the dipole distortions in the distributions of 
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unpolarized quarks in transversely polarized nucleons as seen in impact-parameter space 
and momentum space. In this case the correspondence is between the average distortion 
measured through the anomalous magnetic moment k and the Sivers function. For there 
exist also lattice calculations [74] which refer to a renormalization scale fi"^ = 4 GeV^, and 
give the results = 3.0 and = 1.9. The instant-form version of the xQSM of Ref. [75] 
give = 3.56 and = 1.83. Considering the ratio of up to down contribution, which 
is renormalization scale independent, we find a very good agreement between the lattice 
calculations (1.51) and our model results (1.53 in the LCCQM and 1.48 in the xQSM). 

4 Conclusions 

In this work we presented a first study of GTMDs, which are quark-quark correlators where 
the quark fields are taken at the same light-cone time. By taking specific limits or pro- 
jections of these GTMDs, they yield PDFs, TMDs, GPDs, FFs, and charges, accessible in 
various inclusive, semi-inclusive, exclusive, and elastic scattering processes. The GTMDs 
therefore provide a unified framework to simultaneously model these different observables. 

We took a first step in this modeling, by considering a light-cone wave function (LCWF) 
overlap representation of the GTMDs and by restricting ourselves to the 3Q Fock compo- 
nents in the nucleon LCWF. At twist-two level, we studied the most general transition 
which the active quark light-cone helicity can undergo in a polarized nucleon, correspond- 
ing to the general helicity amplitudes of the quark-nucleon system. We develop a formalism 
which is quite general and can be applied to many quark models as long as the nucleon 
state can be represented in terms of 3Q without mutual interactions. For the radial wave 
function of the quark in the nucleon, we studied two phenomenological successful models 
which include relativistic effects : the light-cone constituent quark model (LCCQM) and 
the chiral quark-soliton model (xQSM). In the LCCQM, the nucleon LCWF was approx- 
imated by its leading Fock component, consisting of free on-shell valence quarks. In the 
xQSM, the quarks are not free but bound by a relativistic chiral mean field, which creates 
on the one hand a shift of the discrete level in the one-quark spectrum, and on the other 
hand also a distortion of the Dirac sea. The latter can be interpreted as arising due to the 
presence of additional quark-antiquark pairs in the nucleon LCWF. Therefore, the xQSM 
has the potential to systematically go beyond a description in terms of the leading 3Q Fock 
component in the nucleon LCWF. Restricting ourselves to the 3Q component in the present 
work, the quarks in both models being either free or interacting with a relativistic mean 
field, allowed us to connect light-cone helicity to canonical spin, and thus relate LCWFs 
with equal-time wave functions. It was shown that the resulting boosts, which connect 
both, in general introduce an angle between the light-cone helicity and the canonical spin 
for a quark with non-zero transverse momentum. 

We firstly obtained the various TMDs in the forward limit, involving either no helicity 
flip, one unit of helicity flip of either quark or nucleon, or transitions where both quark 
and nucleon helicities flip. As the helicity flip processes are accompanied by a change of 
orbital angular momentum, the TMDs reveal such orbital angular momentum in the nu- 
cleon LCWF. We provided predictions for the various TMDs in both models, within the 
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3Q framework, and found several relations between TMDs for both models. The amount 
of orbital angular momentum in the LCCQM in general was found to be larger than in the 
xQSM. In particular, the TMD where both quark and nucleon helicities flip (pretzelosity) 
was found to be about twice as large in the LCCQM as compared with the xQSM. 

We next applied our formalism to describe the GPDs entering hard exclusive processes. 
We found a strong analogy in the multipole patterns in TMDs and GPDs, where the role of 
the quark transverse momentum k± in TMDs is played by the momentum transfer to the 
hadron in the case of GPDs. We found a qualitative difference between both models 
for PDFs or GPDs at x = (for ^ = 0), or in general for GPDs at x = S^, corresponding 
to zero longitudinal momentum for the active quark in the final state. In the LCCQM all 
distributions vanish at this point, in contrast to the xQSM. We also compared PDFs in 
both models, and found that the non-singlet (valence) components of unpolarized (polar- 
ized) PDFs /i (gi) compare reasonably with the phenomenological extractions after NLO 
evolution to 5 GeV^, but systematically undershoot the large-x tail for the tt-quark dis- 
tribution. For the transversity distribution, our results were found to be larger than the 
available parameterizations for both up and down quarks. 

We also compared the results obtained in both models for form factors and charges. 
For the electromagnetic FFs of the proton, and the neutron magnetic FF, the LCCQM was 
found to reproduce relatively well the dependence up to about ~ 1 GeV^, whereas 
the xQSM shows too small magnetic radii. For the neutron electric FF, the xQSM gives 
here a description in good agreement with the data, whereas the LCCQM falls short of the 
data, by about a factor 2 in the range up to ~ 1 GeV^. This behaviour can be improved 
by considering SU{6) breaking terms or higher Fock components, describing the physics of 
the pion cloud. Finally we found for the isovector axial charges as values = 1-24 (1.44) 
within the LCCQM (xQSM) respectively. Furthermore, the quark helicities carry 75 % (85 
%) of the nucleon helicity in the LCCQM (xQSM) respectively. For the tensor charges, 
both LCCQM and the xQSM results evolved at LO to = 0.8 GeV^ agree well for the 
down quark with the available phenomenological extraction, whereas the model result for 
the up quark is larger by ~ 40% compared to the phenomenological value. 

As discussed above, the presented LCWF overlap framework for GTMDs can be sys- 
tematically extended beyond the 3Q LCWF Fock component. Such an extension to the 
5-quark Fock component and its application on the different observables discussed in this 
work, is an interesting subject for future work. 
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A Light-Cone Kinematics 



In this paper, we work in the symmetric infinite momentum frame, where P"*" is large, 
P± = 0± and A • P = 0. The four-momenta involved are then 



P 
k 
A 

n 



-2eP+,2eP",A_ 
0, ±1,0x1 , 



(A.l) 



with P- 



M^+A^ /4 

2(i-gi)'p+ • Note that the form used for n is not the most general one but 
leads to an appropriate definition of TMDs for semi-inclusive deep inelastic scattering and 
Drell-Yan processes. 

For the active quark (i = 1) the momentum coordinates in initial and final hadron 
frames are then given by 



^1 = (y'l 5 



1 - xi Ai 
1-xi A± 



(A.2) 



1 - ^ ' "'^ ' 1 - ^ 2 
while for the spectator quarks {i = 2, 3) they are given by 

Xi A 



ki — {yi, Ki±) 

K = iyi,4±) 



Xi J* 

ki± + 



1 + C 

Xi 



l + C 2 



(A.3) 



1 2 



In Eqs. (A.2) and (A.3), ki± is the average quark transverse momentum and Xj the average 
quark longitudinal momentum fraction in the symmetric frame. 

B LCCQM and xQSM 

We collect in this appendix all the necessary material for the model evaluations. The 
normalizations of the wave functions are fixed by requiring that there are two up and one 
down quarks in the proton. 

B.l LCCQM 

In the numerical evaluations within the LCCQM, we adopted a power-law form for the 
momentum wave function ■0(r) [76] 

i^{r) = 2(27r)' 



1 



yiy2y3Mo {Ml + f3^p' 



(B.l) 
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with the parameters (3 = 0.607 GeV and 7 = 3.5, and the quark mass m = 0.263 GeV. 
These parameters have been fitted to reproduce at best the proton magnetic moment and 
the axial charge. Since 7 > 2, this wave function clearly vanishes as any yi — )• 0, see 




Figure 12. Power-law 3Q LCWF ^{r) used in the LCCQM as function of j/i and y2 with kix — 0.2 GeV 
and K2X = 0.4 GeV. 

figure 12. In eq. (B.l) and in figure 12 the normalization factor is left out. The intrinsic 
scale of the LCCQM is considered to be about Q'q = 0.259 GeV^, see section 3.3. 

B.2 xQSM 

The one-quark discrete-level wave function in the xQSM [21] appears as the sum of a 
bare discrete- level wave function Fj^j(k) and a relativistic contribution F^'^{k) due to the 
distortion of the Dirac sea 

FUk) = F}^:{k) + FZ%k). (B.2) 

The index A refers to the light-cone helicity of the bound quark and a to its canonical 
spin^. 

The bare discrete-level wave function is given by 

(B.3) 

Act 

where At = |k| and the z component of the three- vector k is given by = yAiiy — ^icv with 
the classical soliton mass M-n = 1.207 GeV and the energy of the discrete level E'lcv = 0.2 
GeV. The functions h{K) and j{K) are the upper and lower components, in momentum 

^In Ref. [21] one writes the quark wave function as F-''^ , where j is the quark isospin and a the light-cone 
helicity (merely referred to as "spin" ) . In the model, canonical spin s and isospin j are coupled into grand- 
spin equal to zero. This means that any rotation in spin space can be compensated by a rotation in isospin 
space. One can then interchange canonical spin and isospin indices provided a multiplication by ejs- Our 
quark wave function is then related to the quark wave function of Ref. [21] as Fxa^ = \ ^^^^ ej^F^'^ . 



h{K) + {Kz + iK± X CJ_l) 
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space, of the Dirac spinor for a bound quark. Up to a global normalization factor, we 
found that they are very well approximated by the following forms 

l+AjK + Bj 



h{K) 



A, 



l + BhK^ 

with the parameters 

Ah = 0.64, 
Bh = 12.77, 
Ch = 1.11, 



and 



l + CjK + Dj 



K e 



(B.4) 



0.72, 
4.20, 
-1.06, 



12.40, 

3.48, 

1.17, 



(B.5) 



all in appropriate GeV units (see figure 13 for a comparison with the exact numerical 
solutions). 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 13. Comparison between the exact numerical functions h{hi) and (thin sohd) and their 

approximate form given by eq. (B.4) (thick dashed). The normahzation factor is left free. 

The relativistic contribution to the discrete-level wave function due to the distortion 
of the Dirac sea is given by 



d?k'Wx.,X'a'{kk') 



(J3 h{K ) — K ■ a — — 



(B.6) 



where d^k' = dy' d^K^/(27r)^. It depends linearly on the quark-antiquark pair wave func- 
tion Wxa,x'a'{k,k') whose approximate expression is 



Wxa,x'a'{k,k') 



MqMn 
2ttZ 



Sa<7,av {k, k') - lixa^X'a' {k, k') 



T,xa,x'a'(k, k') = 5^/^ Mqi^y' - y)a-i + Q±-aA 



XX' 



(B.7) 



nAcr,A'cr'(^'^') 



n(g) 



(q-a)^,^ -MQ{y + y') + iQ_Lxa_L 



XX' 



where Mq = 0.345 GeV is the constituent quark mass and q = + {y + y')-M.N) 
is the total momentum of the quark-antiquark pair. In order to simplify the notation we 
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Figure 14. Scalar (dashed blue) and pseudoscalar (solid red) parts of the relativistic mean field in the 
baryon. 

used 

Q± = - y'^± and Z = yy'M at [g^ + w + uj'] (B.8) 

with u = (M^ + K\)/yMN and w' = (M^ + k'D/v'Mn- The functions and n(g), 
shown in figure f4, correspond respectively to the scalar and pseudoscalar parts of the 
relativistic mean field in the baryon. Note that, in principle, the constituent quark mass 
Mq(p) in the model depends on the quark momentum p. This can be seen as a form factor 
that cuts off momenta at some characteristic scale which corresponds in the instanton 
picture to the inverse average size of instantons 1/p ~ 0.6 GeV. One then usually considers 
the scale of the xQSM to be about Qq = 0.36 GeV^. However, since in this study we 
restrict ourselves to the 3Q component only, we may consider that the effective scale of our 
calculations is actually about Qq = 0.259 GeV^ like in the LCCQM. In actual calculations, 
the constituent quark mass is replaced by a constant Mq = Mq (0) and the decrease of the 
function Mgip) is mimicked by the UV Pauli-Villars cutoff at Mpv = 0.557 GeV [77, 78]. 
This value has been chosen from the requirement that the pion decay constant F-,^ = 93 
MeV is reproduced from Mq = 0.345 GeV. The integrals in eq. (B.6) are convergent and 
do not require regularization. Nevertheless, it is not clear to us whether one should still 
apply the Pauli-Villars prescription or not. For this reason, when studying the effects of 
the relativistic contribution F^^'^, we performed calculations both without and with Pauli- 
Villars prescription. 

The discrete-level wave function Fxaik) seems to be a quite complicated function of k, 
A and a. Fortunately, its structure is a bit simpler than at first sight 

FUk)=( ^z/^;;"^ /y^^^'^^^V (B.9) 

where k^^l = Kx iny. The discrete-level wave function has just two independent com- 
ponents f//{y,K±) and f±{y,K±) which are functions of y and k± only. If we neglect the 
relativistic contribution F^^{k), these functions can be written as 

f//{y,K±) = h{K) + Kzj{K)/K and f±{y, k±) = (B.IO) 
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Three-dimensional plots of these functions are shown in figure 15. Unless mentioned ex- 
plicitly, all the results presented in this paper for the xQSM have been obtained with this 
wave function. One expects the relativistic contribution F^'^{k) to change Eq. (B.IO) in 




Figure 15. The two independent components f//{y,Kj_) and f±{y,ii±) of the bare discrete-level wave 
function They are normalized so as to have two up and one down quarks in the proton. As long as we 

neglect the relativistic contribution F^''^{k), f±{y, k±) is actually a function of k = ^y (yM m — -Eiov)^ + n\ 
only, as one can see from Eq. (B.IO). The thick elliptic curves represent constant values of k. 

an appreciable manner only around y = 0. This is confirmed by an explicit calculation 
of the complete discrete- level wave function Fx„{k), both without and with Pauli-Villars 
prescription. In figure 16 the discrete-level wave function in three versions of xQSM are 
compared. 

B.3 Generalized Melosh rotation 

There are striking similitudes at the 3Q level between the xQSM and the LCCQM. For ex- 
ample, both models exhibit the SU (6) spin-flavor symmetry and are based on a completely 
symmetric momentum wave function. In the LCCQM the quarks are free while they are 
bound by a relativistic mean field in the xQSM. This means that in both cases the quarks 
are considered without mutual interactions. Usually, it is an incredibly difficult task to 
relate the (instant-form) canonical spin to the light-cone helicity. But since in the models 
we considered here there are no mutual interactions, the relation turns out to be rather 
simple. In the LCCQM, the quarks being free, it is well known that their canonical spin 
and light-cone helicity are just related by a Melosh rotation [25], which is a special case of 
Wigner rotation. The amplitude and the axis of rotation depend on the quark momentum. 
Is it the same in the xQSM? Since here the quarks are not free, it is clear that one cannot 
use the Melosh rotation. In this model, it is the discrete-level wave function F\„{k) that 
tells us how to relate the canonical spin a to the light-cone helicity A, and that this relation 
depends on the quark momentum k. Dividing Fxa-{k) by f'j + i'^'±f± in Eq. (B.9), we obtain 
a unitary 2x2 matrix, i.e. here also light-cone helicity and canonical spin are just related 
by a rotation. This explains the introduction of the generalized Melosh rotation given by 
Eq. (2.17). 
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Figure 16. Comparison between different versions of the discrete-level wave function, see text. Solid red 
curves: only (xQSM I). Dashed blue curves: + F""'' (xQSM II). Dotted black curves: F'°^ + F'""^ 
with Pauli-Villars regularization (xQSM III). The light bands give the range of possible values for the 
discrete level, depending on the treatment of the sea contribution. Clearly, only the region around y = is 
affected in an appreciable manner by the relativistic contribution F'^^^ , 

If we consider the limit of vanishing mean field in the xQSM, the quarks become free 
and we expect the generalized Melosh rotation to reduce to the ordinary Melosh rotation. 
Let us check this explicitly. The relativistic contribution F^"^'"^ is directly proportional to the 
mean field represented by and n(q'). This means that in the limit of vanishing mean 
field F^^'"^ — )■ 0. Let us then focus on Since quarks becomes free, we have to replace 

in Eq. (B.3) the discrete-level energy -Bicv) the soliton mass M.n and the constituent quark 
mass Mq by the free quark energy w, the free invariant mass A^o and the free constituent 
quark mass m, respectively. Moreover one obtains that the upper and lower components 
are simply proprotional j{K) — )• Kh{K)/{uj + m). Collecting all the parts, we find 

FxAk) ^ Hk) [m + yMo + ia-{e, x k^)]^^ = f{k)D'J^*{R,f{k)), (B.ll) 

where f(k) is some momentum wave function and D^J^* {Rcfik)) is the matrix element 
of the Melosh rotation, see Eq. (2.18). As announced, we naturally recover the Melosh 
rotation in the limit of vanishing mean field. 

C LCWF overlap 

We derive in this appendix the master formula written in section 2.5. Using the expression 
for the LCWF given by Eq. (2.16) for the overlap representation of the correlator tensor 
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W'^'^ , we obtain Eq. (2.20) with the tensor A'^'^ given by 
■^'•"{r'^r) = 1 E i^n"-^' ^t''^' ^T'''' K'^.S^'iM^ 



(C.l) 



where N'^{k',k) = D'^ {k') D{k) are 2x2 matrices. We decompose D{k) and N^{k',k) 
on the basis = and introduce the corresponding components and O^'^ 



D{k) = a^d^', df'=(e,-K,ie,xK 
N^{k\k) = ^T^O^^ 



(C.2) 
(C.3) 



where we used the notation a = d/\d\. Since |Tr [cj'^fj'^] = g^'^ , the components of the 
matrices N"{k' , k) in the basis are given by 
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= iT,- 

2 



a>'D^k')a''DCk) 
d'^dp. 



„(y- -V „B 



(C.4) 



Now using the identity (eoi23 = +1) 



It,- 

2 



^tJ- T^Ol 



we obtain 



O^'" = d'^^'d" + - [d'* ■ d) g^" + ie^"''"^d'*d^ 



or more exphcitly^ 
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\K'\\K\ 



K'-K i K' X K] i K'xK 



i K' X K 



(C.5) 
(C.6) 



i{K' X K) K' -K - 2K'^K^ -K^y - K'yK,^ K!^K^ + K'^K, 



i[K' xK] -K'K^ - K'^Ky K'-K- 2K'Ky K'K, + K'.Ky 



yJ.^y ±^y2 

i(k' X K) -K'^K^ - K'^K^ -KKy - K'yK^ -K' ■ K + 2K'^K^ 



Let us come back to A^'^ . We can now write 



(C.7) 



A',A,(7',o-, 



*Note that the elements of this matrix have already been obtained some years ago using the Melosh 
rotation in the LCCQM. With the notations of refs. [13-15], the matrix C^" reads 
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where we defined the four-component vectors = . We have to sum over all the 
polarizations A',A, C7^,c7i which seems a priori quite involved. We can however already 
guess the actual form of the result by standard tensorial analysis. Indeed we must end up 
with a tensor with indices and v, and this tensor has to be constructed out of Of^, I2 
and Z3. There are only three possible terms and A^^ then takes the form 

A^'^ = AO^^" {h ■k) + B It^ {k ■ OiY + Cl^ {h ■ OiY . (C.9) 

Since spectator quarks are equivalent, we should have B = C . If one is not interested in 
the quark flavor, then all the three quarks are equivalent and we have A = B = C which 
can be absorbed in the normalization of the LCWF. 

For the purpose of the present paper, let us determine the flavor coefficients A and B 
in the case of a proton target. Since in the process considered each individual quark flavor 
is conserved, it is convenient to divide the SU{Q) spin-flavor wave function (J)^i°"20'3 -j^^^ 
three terms 

^ K,uud ^ ^ k,udu ^ ^A,duu ■ \^\-j.±\jj 

Up to an overall normalization factor, ^^^^^^^ can be written in terms of Kronecker and 
Levi-Civita symbols (see e.g. [21, 23]) 

(f,o-io-2(T3 _ xo-i 0-20-3 I ro2 0103 (n ^^\ 

The udu and duu terms are obtained from the uud term by permuting the labels (2, 3) 
and (1,3), respectively. Using the relation (a^)^^'^ e^'^' e"^^ = (a^)^^' and then the iden- 
tity (C.5), we get the coefficients 





AP 


BP 


CP 


uud 
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-1 


udu 
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-1 
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duu 


-1 
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2 



(C.12) 



Since the first quark was chosen to be the active one, the sum of the two first lines corre- 
sponds to the contribution of the up flavor, while the last line corresponds to the contri- 
bution of the down flavor. 
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